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I 1 ABSTRACT 

Q Context. The physical interpretation of scattering line polarization offers a novel diagnostic window for 
^ exploring the thermal and magnetic structure of the quiet regions of the solar atmosphere. 

Aims. To evaluate the impact of isotropic colUsions with neutral hydrogen atoms on the scattering polariza- 
tion signals of the 13 lines of multiplet 42 of Tii and on those of the K Une and of the IR triplet of Can, 

T— I 

^ with emphasis on the colUsional transfer rates between nearby 7-levels. 

Methods. We calculate the Unear polarization produced by scattering processes in a plane-parallel layer 
illuminated by the radiation field from the underlying solar photosphere. We consider realistic multilevel 
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models and solve the statistical equilibrium equations for the multipolar components of the atomic density 
matrix. 

Results. We give suitable formulae for calculating the collisional rates as a function of temperature and hy- 
drogen number density. We confirm that the lower levels of the 1 3 lines of multiplet 42 of Ti i are completely 
depolarized by elastic colUsions. Upper-level collisional depolarization turns out to have an unnoticeable 
^ impact on the emergent Unear polarization amplitudes, except for the /14536 line for which it is easier to 
notice a rather small depolarization caused by the colUsional transfer rates. Concerning the Can Unes, we 
show that the collisional rates play no role on the polarization of the upper level of the K Une, while they 
have a rather small depolarizing effect on the atomic polarization of the metastable lower levels of the Can 
IR triplet. 

Conclusions. Although the colUsional transfer rates between nearby 7-levels seem to play a minor role for 

most of the lines we have considered in this paper, except for the magnetically insensitive /14536 line of 
Til, they might be important for other atomic or molecular systems with closer y-levels (e.g., hyperfined 
structured multiplets and/or molecules). Therefore, future research in this direction will be worthwhile. 



Key words. Scattering - Polarization - Atomic processes - Sun: photosphere - Sun: chromosphere - Line: 
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1. Introduction 

The physical interpretation of observations of the spectral line polarization produced by scatter- 
ing processes offers a unique opportunity to learn about the thermal and magnetic properties of 
the "quiet" regions of the solar atmosphere (e.g., the recent reviews by Stenflo 2003 and Trujillo 
Bueno 2003). Of great scientific interest is the possibility of obtaining empirical information on 
hidden aspects of solar magnetism, such as tangled magnetic fields in the solar photosphere (e.g., 
Trujillo Bueno et al. 2004; Stenflo 2004; Manso Sainz et al. 2004), and/or very weak magnetic 
fields in the solar chromosphere (e.g., Manso Sainz & Trujillo Bueno 2001, 2003, 2007; Trujillo 
Bueno & Manso Sainz 2002; Holzreuter et al. 2006). In order to be able to achieve such goals in 
a reliable way it is crucial to carry out detailed radiative transfer calculations using sufficiently 
accurate values of the rates of isotropic collisions with neutral hydrogen atoms, such as those 
computed by Derouich et al. (2003; 2004; 2005a, b) within the framework of the semi-classical 
theory of collisions. The present paper addresses this issue by investigating their influence on the 
scattering polarization of some spectral lines of the linearly polarized solar limb spectrum, which 
are of great diagnostic value. 

Collisions with neutral hydrogen atoms in the solar atmosphere may modify the atomic level 
polarization that anisotropic radiation pumping processes induce in the atomic and molecular 
energy levels. We have elastic collisional rates, D^^\j), which normally tend to destroy the 
population imbalances between the magnetic sublevels pertaining to each level of total angular 
momentum J. In addition, we may also have significant transfer rates between diff'erent J levels 
due to inelastic and super-elastic collisions, C^f^\j, Ji) and C^f \j, Ju), respectively. These col- 
lisional transfer rates may modify the population and/or the atomic polarization of the level J 
under consideration through colhsional transfer from the lower (7/) and upper (/„) levels. Here 
we pay particular attention to both type of collisions with neuttal hydrogen atoms, with empha- 
sis on the (typically neglected) collisional transfer rates between the / levels pertaining to each 
spectral term of the chosen multilevel model. 

In this paper we focus our attention on the 13 lines of multiplet No. 42 of Ti i and on the 
following lines of Ca n: flie K Une at 3934 A and flie IR triplet at 8498 A, 8542 A and 8662 
A. The reason for choosing multiplet 42 of Ti i, which results from the ttansitions between the 
upper tsTmy^F^ and the (metastable) lower term a^F, is because one of the lines (i.e., /14536) is 
insensitive to magnetic fields (the Lande factors of its lower and upper levels are zero), while the 
remaining 12 lines are sensitive to the Hanle effect. Therefore, ratios of the observed polarization 
amplitudes in suitably chosen line pairs can be used to infer the presence of tangled magnetic 
fields in the quiet solar atmosphere (Manso Sainz et al. 2004). The scientific motivation for 
choosing the above-mentioned lines of Ca ii is because they are of great diagnostic interest for 
investigating the thermal and magnetic structure of the quiet solar chromosphere (Manso Sainz & 
TrujiUo Bueno 2001, 2003, 2007; TrujiUo Bueno & Manso Sainz 2002; Holzreuter et al. 2006). 

* Consejo Superior de Investigaciones Cientfficas (Spain) 
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2. Formulation of the problem 

We consider a plane-parallel layer, of Ti i atoms or Ca ii ions in the absence of magnetic fields, 
illuminated by the radiation field of the underlying solar photosphere. We assume that this 
anisotropic radiation field is axisymmetric with respect to the local vertical and unpolarized. 
Therefore, the following spherical tensors components suffice to describe the optical pumping 
properties of the incident radiation field (e.g., Landi Degl'Innocenti & Landolfi 2004) 
1 

J« = - J ^ mdM, (1) 



-1 

jl = — \ [(3i/ - i)mw, (2) 

4 V2 J-i 



where is the intensity and /u the cosine of the angle that the radiation beam under considera- 
tion forms with the local vertical. In the reference system where the z-axis (i.e., the quantization 
axis of total angular momentum) is chosen along the symmetry axis of the incident radiation field 
(i.e., along the solar local vertical direction), the anisotropic illumination of the atoms produce 
population imbalances which we quantify by the following multipolar components of the atomic 
density matrix: 



M=-J 



J J 
M -M 



pj(M,M), (3) 



where pj{M, M) is the population of the substate with magnetic quantum number M, while K = 
0, 2, 2J for the Ti i levels (Ti i levels have integer total angular momentum values) and K = 
0, 2, 27-1 for the Ca n levels (Ca ii levels have half-integer total angular momentum values). 
Note that Pq(/) is proportional to the total population of the J level under consideration, while 
p^{J) is the so-called alignment of the level (which quantifies the degree of population imbalances 
among its magnetic sublevels). Note also that we are safely neglecting the hyperfine structure of 
the odd isotopes, given that the most abundant isotopes have zero nuclear spin (87% for titanium 
and practically 100% for calcium). 

Both, radiative and coUisional transitions contribute to the Pq (/) values corresponding to 
each J level of the assumed atomic model. In the regime of the impact approximation the rate of 
change of Pq{J) can be written as 

J^PqU) = [|po(-/)]rad + [j^P^(J)]coU. (4) 

The contribution of the radiative rates has been studied in detail elsewhere (e.g., Manso Sainz 
& Landi Degl'Innocenti 2002; Manso Sainz et al. 2006) and its importance will not be empha- 
sized further here. Adopting the same notation of the monograph by Landi Degl'Innocenti & 
Landolfi (2004) we write the contribution of collisional transitions as follows: 
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+ ^ Cf\a„J,„ aJ) + D^^XaJ)] p^{aJ) (5) 

aiJi 
a I, J II 

where D'^^HaJ) (with D^^HflJ) - 0) are the rates of elastic collisions that induce transitions 
between magnetic sublevels belonging to the same \aJ > level, C\^\aJ, aiJi) denote the multi- 
polar components of the inelastic coUisional rates that induce transitions between a lower level 
\aiJiMi > and the level \aJM >, and C^^\aJ, auJu) are the multipolar components of the supere- 
lastic collisional rates that induce transitions between the upper level |a„y„M„ > and the level 
\aJM >. 

Given that in this paper we are especially interested in the collisional transfer rates between 
the 7-levels of each individual term, we actually have that a - ai - a„ in Eq. (5). CoUisional 
transitions between the lower and upper levels of the spectral Unes under consideration induced 
by inelastic collisions with electrons are disregarded here (i.e., we assume that the atomic excita- 
tion is dominated by radiative transitions). Therefore, here we are investigating a problem where 
the radiative transitions between the lower and upper levels of the spectral lines under consid- 
eration induce population imbalances among the magnetic sublevels of such levels, which are 
modified by the presence of isotropic colUsions with neutral hydrogen atoms, both by the elastic 
colUsions and by the collisional transfer rates between "nearby" 7-levels only (i.e., between the 
7-levels belonging to each spectral term of the multilevel model under consideration). 

Note that a different notation is used in the previous papers by Derouich et al. (e.g., Derouich 
et al. l2004l l. There the symbol D^iaJ — > a J') is used to denote the population transfer rates from 
J to 7', while D^{aJ — > aJ') refers to the polarization transfer rates. The relationship between 
the two notations is: 

D^iaJ^aJ')^ ^J^^ C'^'^\aJ',aJ). (6) 

We obtain the Pq(J) elements corresponding to each /-level of the atomic model under con- 
sideration by assuming statistical equilibrium (i.e., |^ M"-^) J ^ ^ [ ^"j"'''' ] ^ - 0)- Since the 
ensuing algebraic system of equations in the unknowns Pq(J) is singular, one of the equations is 
substituted with that expressing the conservation of particles (i.e., with YjJ P^iJ) V27 + 1 = 1). 
We have caiTied out this type of calculations imposing the above-mentioned anisotropic illumina- 
tion conditions for the following cases: (a) neglecting collisions, (b) considering only the elastic 
colUsions and (c) including also the inelastic and superelastic colUsional transfer rates between 
nearby J levels (i.e., between those pertaining to the same spectral term). 

From the computed p^iJ) values we need to use only the K - and K - 2 components 
corresponding to the upper and lower levels of the chosen spectral line, since these are the density 
matrix elements needed to compute the / and Q components of the emission vector and of the 
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propagation matrix that appear in the Stokes-vector radiative transfer equation (see Eqs. 7.21 in 
Landi Degl'Innocenti & Landolfi 2004). The emissivity contributions, which are proportional to 
the density matrix elements of the upper level, are given by 

= eo[p°o(Ju) + wf\^(3i/ - 1)pI{Ju)1 (7) 

= eowf\ ^ ( 1 - l^^)pliJu), (8) 

where eo - {hvlAn) A^e V2/« + 1 A'^ (N being the number density of the atoms under consider- 
ation), and wf-'j^ is a numerical factor that depends only on the quantum numbers of the levels 
involved in the transition (e.g., see Table 10.1 in Landi Degl'Innocenti & Landolfi 2004). The 
absorption and dichroism contributions, which are proportional to the density matrix elements of 
the lower level, are given by 

= m{Po(J() + ^Tj -^(V - i)pl(Je)l (9) 
' "2 V2 

= ^owfX:^i^ - t^^)pliJeX (10) 

where 770 = (hv/A-n) B(u ^2Ji + 1 N. Note that in the expressions for eg and rjQ the reference 
direction for Stokes Q has been chosen as the perpendicular to the radial direction through the 
observed point. 

As in Manso Sainz & Landi Degl'Innocenti (2002), we calculate the emergent fractional 
linear polarization in the limit of tangential observation in a plane-parallel atmosphere (90° scat- 
tering geometry). For such a hne-of-sight the solution of the transfer equations for Stokes / and 
gis 

iQ\ ^ ^qI^i - riQlm ^^^^ 

I / /;i~.0 1 - (r]Q/r]i)(£Q/€i) ' 

where ex and j]x (with X - I oiX = Q) must be evaluated at the surface point of the plane-parallel 
slab. Actually, given that the solar atmosphere is a weakly anisotropic medium, this expression 
can be simplified as follows: 



X eg/ei - i]q/j]j, (12) 





which shows clearly that the observed scattering polarization in a given spectral line has, in gen- 
eral, two contributions (TrujiUo Bueno 1999): one due to the selective emission of polarization 
components resulting from the polarization of the upper level, and an extra one caused by dichro- 
ism -that is, by the selective absorption of polarization components resulting from the lower level 
polarization. 

Obviously, when the lower level of the spectral line under consideration is unpolarized, then 
the only contribution to the emergent fractional linear polarization would be that caused by the 
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selective emission of polarization components that result from the population imbalances of the 
upper-level. Under such circumstances 



where the approximate expression applies to the line center of a sufficiently strong spectral line. 
In Fig. 2 below we will show what happens when we make this unpolarized lower-level assump- 
tion. In addition, we will also show what happens when the lower level is allowed to be polarized. 
To this end, we only need to solve the statistical equilibrium equations without introducing any 
colUsional depolarization. An important point to remember here is that the presence of lower 
level atomic polarization may have a significant feedback on the polarization of the upper levels, 
so that a calculation of the emergent Qjl neglecting the dichroism contribution (i.e., ignoring the 
tiqItji term of Eq. 12) will give in general a Qjl amplitude different to that corresponding to the 
unpolarized lower level case (Trujillo Bueno 1999, 2001). 

3. Application to multiplet 42 of Ti i 

This section describes the results of our investigation of the influence of elastic, inelastic and 
superelastic collisions on the atomic polarization of the lower and upper levels of the 13 lines of 
multiplet 42 of Ti i, and on the amplitudes of the ensuing emergent fractional linear polarization. 
We begin by considering the sensitivity of the calculated fractional linear polarization to the 
complexity of the assumed multilevel model, justifying why we end up using a 10 level model 
with the five J levels of the lower term and the five J levels of the upper term of multiplet 
42 of Ti I. We then investigate the sensitivity of the emergent linear polarization to the above- 
mentioned collisions. Appendix A gives polynomial approximants for calculating the coflisional 
rates for any combination of hydrogen number density and temperature of the solar photospheric 
plasma. Such collisional rates, which are given here in s result from the application of the 
semi-classical theory used by Derouich et al. (2005b). 

We have considered two atomic models. The most realistic one (hereafter the full model) 
is similar to that used by Manso Sainz & Landi DeglTnnocenti (2002), which has 395 J levels 
and afl the allowed transitions between them. The total number of Pq(J) elements necessary to 
describe the excitation state of the full atomic model is 1564, whose values can be found by solv- 
ing the above-mentioned statistical equiUbrium equations for given illumination conditions (i.e., 
using the 7g and values computed from the center-to-limb variation of the solar continuum 
radiation tabulated by Cox (2000) for each allowed radiative transition|^ In any case, given the 
low degree of radiation anisotropy in the solar atmos phere, one can safely solve the statistical 
equilibrium equations neglecting the Pq{J) elements with K > 2 (e.g., Manso Sainz & Landi 

' Figure 2 of Manso Sainz & Landi Degl'Innocenti (2002) shows the variation with wavelength of the 
anisotropy factor, w = y/lJ^/J^, and of the number of photons per mode (i.e., h = j'^icr I2hv^)) of the solar 
continuum radiation. 




(13) 
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Fig. 1. Grotrian diagram of multiplet 42 of Ti i showing the angular momentum values of the lev- 
els of the lower term a^F and of those of the upper term j'^F". The wavelengths (in Angstroms 
units) of the 13 hnes of this multiplet are also indicated. 

Degl'Innocenti 2002), thus reducing the number of Pq{J) unknowns to 771. For this reason, in 
the appendices we have given the colhsional rates only for = and K = 2. 

The second multilevel model we have used is that shown in Fig. 1, which contains all the 
lower and upper J levels of multiplet 42 of Ti i. The number of Pq (/) elements necessary to de- 
scribe the excitation state of this 10 level model atom is 20 when those with > 2 are neglected. 
It is important to note that the lower levels of the 13 hnes of this multiplet are metastable, so that 
for a given coUisional rate the population imbalances of such long-lived levels are much more 
vulnerable to collisions than the short-lived upper levels of the spectral hnes under consideration. 

Figure 2 shows the results of calculations with the two multilevel models mentioned above, 
neglecting the colhsional rates (i.e., taking into account only the radiative transitions). The up- 
per panel neglects the dichroism contribution caused by the selective absorption of polarization 
components (i.e., the emergent polarization amplitudes were obtained by using Eq. 13, which 
neglects the /ye///, term of Eq. 12). The lower panel was obtained using instead Eq. (13), so that 
it takes into account that selective absorption contribution, in addition to the usual one caused by 
the selective emission of polarization components resulting from the population imbalances of 
the upper level. As seen in Fig. 2 the 10 level model atom significantly overestimates the frac- 
tional linear polarization of the /14536 fine, when the role of collisions is neglected. As pointed 
out by Manso Sainz et al. (2004) this Ti i line is insensitive to the presence of magnetic fields 
b ecause the Lande factors of its upper and lower levels are exactly zero. Therefore, it is a very 
useful reference line whose polarization pattern can be used to investigate the rehabihty of the 
thermal and dynamic properties of solar photospheric models (see Shchukina & Trujillo Bueno 
2007). 
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Fig. 2. Emergent fractional linear polarization amplitudes, QII, for each of the 13 lines of mul- 
tiple! 42 of Ti I indicated in Fig. 1. □: 10 levels model, x: full multilevel model. Collisions are 
not taken into account. Panel (a) shows the emergent fractional polarizations if dichroism is ne- 
glected (i.e., using Eq. 13). Panel (b) shows the emergent fractional polarizations if dichroism is 
taken into account, in addition to the standard contribution of selective emission of polarization 
components (i.e., using Eq. 12). Note that in panel (b) the polarization amplitudes corresponding 
to the calculation performed using the full multilevel model agree with those shown by Manso 
Sainz et al. (2004) in their Fig. 2b. 

As shown in Fig. 3, when the lower metastable levels are assumed to be completely un- 
polarized then the 10 level model atom gives excellent results. Interestingly, as pointed out 
by Shchukina &. Trujillo Bueno (2007), this seems to be indeed the case in the solar atmo- 
sphere. In fact, the elastic collisional rates given in Appendix A for the 5 lower levels of 
multiplet 42 of Ti i are such that 6(Ji) - D^^\Ji) fiife(//) » 1 for each lower-level Ji, even 
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Fig. 3. Emergent fractional linear polarization amplitudes, Q/I, for each of the 13 lines of mul- 
tiple! 42 of Ti I indicated in Fig. 1. □: 10 levels model, x: full atomic model. The atomic polar- 
ization of all (long-lived) lower levels is assumed to be zero. The influence of collisions on the 
population imbalances of the upper levels is not taken into account. 



when using the relatively low hydrogen number densities of the temperature minimum re- 
gion of standard semi-empirical atmospheric models (i.e., n/^ a; 2x10'^ cm"-'). For example, 
l/fiife(^/) ~ B/«^o ~ BiuBy{T)l2 ~ 1.7x10^ s"' (when using A = 4536A and T = 5800 K). On die 
other hand, D^^\Ji) ~ 8x10"'' «h (see Appendix A). Therefore, 6(Ji) ~ 1 for ~ 2x10'^ cm"\ 
while 6{Ji) ~ 10^ in the denser solar atmospheric region where the line-center polarization of 
such Ti I lines originates. Accordingly, from now on we will use the above-mentioned 10 level 
model atom without lower-level polarization in order to investigate the impact of collisions on 
the population imbalances of the 5 upper levels of Fig. 1, and on the emergent linear polarization 
in the ensuing 13 lines. Finally, it is also of interest to note the differences between the polar- 
ization amplitudes shown in Figs. 2a, 2b and 3. In order to understand why they are different, it 
suffices to take into account what we have mentioned in the last paragraph of Section 2. 

In order to investigate the impact of the collisional rates on the emergent fractional linear 
polarization amplitudes we have done the following. For each spectral line we use the hydrogen 
number density and the kinetic temperature values corresponding to the height h in the quiet Sun 
photospheric reference model of Maltby et al. (1986) where the line-center optical depth is unity 
for a simulated observation along a line-of-sight specified by = 0.1 (with // the cosine of the 
heliocentric angle). Such heights h for each line of multiplet 42 of Ti i are given in Fig. 5 of 
Trujillo Bueno et al. (2006), where it is shown that they are contained in the range between 400 
and 550 km. For example, hx45Q km for the /14536 line (Shchukina & Trujillo Bueno 2007). 
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Figure 4 shows the calculated fractional polarization ampUtudes in the presence of colhsions. 
Elastic colhsions with neutral hydrogen atoms (i.e., the D'^^\j) rates) have a negligible influence 
on the upper levels of aU the lines of the titanium multiplet under consideration, while the col- 
lisional transfer rates have a smaU but noticeable depolarizing impact on a few of them (i.e., 
mainly on the upper levels with = 1 and = 2). This is because the population transfer 
rate between the = 1 and = 2 levels is substantiaUy larger than the rest of the coUisional 
transfer rates (see Appendix A.4). As seen in Fig. 4, the observable efl'ect of such an atomic level 
depolarization is mainly seen in the /14536 hne (i.e., that of Fig. 1 with = = 1), in spite of 
the fact that the /14544.7 line (i.e., that with = 1 and Ji = 2) has the same upper level. For both 
lines we find 



(6/^)all coUisional rates = 75 (14) 

(e/')only the colhsional elastic rates 

but in Fig. 4 this is only noticed in the 2// of the /14536 line. The reason for this behaviour can be 
easily understood through Eq. (13), which implies that A(2//)«(3/2>^)w^^|, A[p^(7„)/pO(7^)], 
with wfj' = -0.5 for the /14536 hne and wf^ = 0.1 for the /14544.7 hne. Our conclusion is that 
the coUisional transfer rates among the upper / levels of multiplet 42 of Ti i have a small but 
noticeable depolarizing impact on the emergent fractional linear polarization of the /14536 hne. 
As mentioned above, this is a truly remarkable spectral hne because its polarization turns out to 
be virtuaUy insensitive to the presence of magnetic fields (see Fig. 2 of Manso Sainz et al. 2004). 
We have checked that the aligrmient transfer rates among the upper levels of the multiplet do not 
change this conclusion (see, e.g.. Fig. 4). 

4. Application to Ca ii 

This section describes the results of our investigation of the influence of colhsions on the atomic 
polarization of the levels of the Ca n model of Fig. 5. In addition to the elastic collisions, we also 
take into account the coUisional transfer rates between the two J levels of the term (which are 
the upper levels of the H and K fines) and between those belonging to the term (which are 
the lower levels of the Ca n IR triplet). In contrast with the lower levels of multiplet 42 of Ti i 
the two metastable levels of the Ca n IR triplet (i.e., the levels with 7 = 3/2 and J = 5/2) 
are polarized. This was demonstrated by Manso Sainz & Trujillo Bueno (2003) who solved the 
fuU radiative transfer problem for the multUevel model atom of Fig. 5 in a semi-empirical model 
of the solar atmosphere, and taking into account elastic, inelastic and superelastic colhsions with 
the ensuing rates calculated foUowing Shine & Linsky (1974) for the inelastic and superelastic 
rates and Lamb & ter Haar (1971) for the elastic rates. Appendix B gives the coUisional rates 
obtained by Derouich et al. (2004), but using instead the notation adopted in the present paper. 
Such rates were obtained via the appfication of the semi-classical theory of collisions. 

In order to investigate the impact of colhsions on the Ca n lines we have appUed the foUowing 
strategy, which relies on the results that Manso Sainz & Trujillo Bueno (2003) obtained through 
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Fig. 4. Emergent fractional linear polarization amplitudes, Qjl, for each of the 13 lines of multi- 
plet 42 of Ti i using the 10 levels model atom shown in Fig. 1 without lower-level polarization, 
x: Qjl amplitudes taking into account only D^^\J). □: Q/I values using all collisional rates: 
C^°\aJ',aJ) and C^^>(aJ,aJ'). 

Table 1. Some of the values used in the Ca ii calculations 



line 


H 


K 


8498 


8542 


8662 


h (Km) 


2040 


2000 


1200 


1460 


1380 


Hh (cm"') 


14 X 10" 


1.7 10" 


7.6 X 10'- 


1.7 X 10'2 


2.5 X 10'^ 


r(K) 


8780 


8550 


6475 


6930 


6800 


n 


0.12 X 10"' 


0.13 X 10"' 


13 X 10"' 


8.7 X 10"' 


9.1 X 10"' 


w 


0.031 


0.031 


0.012 


0.028 


0.036 



detailed multilevel radiative transfer calculations in a semi-empirical atmospheric model similar 
to that of Fontenla et al. (1993). For each spectral line indicated in Fig. 5 Table 1 gives the total 
hydrogen number density, the kinetic temperature, n (i.e., the number of photons per mode of 
the radiation field) and w - V2/q//q (i.e., the anisotropy factor) at the atmospheric height in the 
plane-parallel semi-empirical model where the line center optical depth is unity for a simulated 
observation at /i = 0.1. We performed four multilevel calculations in order to compute the emer- 
gent Q/I in each of the four spectral lines that can in principle be polarizecj^ For example, for 
the AS542 line we solved the statistical equilibrium equations for the 5-level model atom of Fig. 
5 using the rin, T, h and w values corresponding to the height h - 1460 Km. These values at 
h = 1460 Km are given in Table 1 for the /18542 line, but we have used also the n and w values 

^ Note that there can be no scattering polarization in the H line because its lower and upper levels cannot 
be aligned. 
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4s 

1/2 

Fig. 5. Grotrian diagram of Can showing the upper and lower 7-levels of the H and K lines and 
of the IR triplet. 

corresponding to each of the remaining line transitions at the very same height. The resulting pg 
values corresponding to the lower and upper levels of the /18542 hne were then used to compute 
QII according to Eq. (12). A similar calculation was performed for each of the other three lines 
that can be polarized by scattering processes (i.e., the K hne, the 8498 A line or the 8662 A line), 
but solving the statistical equihbrium equations using for each of the hne transitions of Fig. 5 the 
Hh, T, n and w values corresponding to the atmospheric height h where the hne center optical 
depth along ju = 0.1 is unity for the particular spectral hne under consideration. 

The results are shown in Fig. 6 for the coUisionless case, for the case in which only elastic 
colhsions are taken into account, and for the most general situation in which also the coUisional 
transfer rates are considered. As seen in the figure, the coUisional transfer rates play only a minor 
role on the scattering polarization of the Ca n IR triplet, while they are completely insignificant 
for the K line at 3934 A. 

Finally, we would hke to coimnent briefly on an interesting point about the K-line polariza- 
tion. It is commonly beheved that isotropic coUisions have always a depolarizing role. However, 
when the lower levels of the line transitions under consideration are polarized, we may find cir- 
cumstances where isotropic coUisions increase the emergent polarization amphtudes. This was 
first pointed out by TrujUlo Bueno & Landi Degl'Iimocenti (1997) for the two-level atom case 
with lower-level polarization, but other authors have supported this conclusion through their 
calculations for more reahstic model atoms (e.g., Manso Sainz & Landi Degl'Iimocenti 2002; 
Klement & Stenflo 2003). 
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Fig. 6. Emergent fractional linear polarization amplitudes, QII, for the K line and for the IR 
triplet of Ca ii. The calculations have been performed using the multilevel model of Fig. 5, and 
following the strategy indicated in the text. ■: Qjl values in the absence of collisions, x: Qjl 
amplitudes taking into account only D^^\J). □; Q/I values calculated using aU coUisional rates: 

D^^\J), C'-^Kaf, aJ) and C^^\aJ, aJ'). 

Here we show that there is a range of values of the hydrogen number density for which we 
have coUisional enhancement of the emergent fractional linear polarization in the K line of Ca ii. 
As seen in Fig. 7, the hh values for which this is possible are larger than those encountered in the 
region of the solar chromosphere where the K line-core polarization originates, so that collisions 
with neutral hydrogen atoms can indeed be neglected when modeling the K-line polarization. In 
any case, it is important to keep in mind that this possibility exists for some spectral lines of the 
solar spectrum. 

5. Concluding comments 

The coUisional rates given in the appendices resulted from the application of the semi-classical 
theory of collisions. From comparisons with the published results obtained through quantum 
chemistry methods we estimate that any possible error in our colUsional rates for Ca ii is smaller 
than 10%. Unfortunately, we cannot give a similar error estimate for our coUisional rates of the 
complex Ti i atom. The only thing we can say at this stage is that the semi-classical theory of 
colUsions appUed to levels where we have just one single valence electron is expected to give 
fairly accurate results for the intermediate interatomic distances, which are the ones playing the 
key role in the calculation of the coUisional rates (see Fig. 6 of Derouich et al. 2003). 

We confirm that elastic colUsions with neutral hydrogen atoms in the solar atmosphere de- 
stroy the atomic polarization of the (metastable) lower levels of multiplet 42 of Ti i. This is im- 
portant because without lower- level polarization a 10-level model atom is sufficient for modeling 
the scattering polarization of the 13 lines of such a multiplet via radiative transfer calculations 
in three-dimensional (3D) models of the solar atmosphere, as done by Shchukina & Trujillo 
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Fig. 7. Emergent fractional linear polarization amplitude of the Ca ii K line as a function of the 
neutral hydrogen density, (a) Dot-dashed line: Q/ 1 values calculated assuming that the metastable 
levels ^Z)3/2 and ^1)5/2 are completely unpolarized. (b) Solid line: Q/I amplitudes calculated us- 
ing all collisional rates, taking into account the polarization ofsuch metastable levels: D^^\J), 
C^"\aJ',aJ) and C^^\aJ,aJ'). (c) Dotted line: same as in (b) but taking into account only 
D^^\j). The dashed vertical lines, between hh - 2 x 10'^ cm""* and hh - 2 x 10'"* cm""*, 
represent the collisional regime where the scattering polarization in the K line increases with the 
hydrogen number density. 



Bueno (2007). On the other hand, elastic collisions have a rather small depolarizing effect on the 
(metastable) lower levels of the Ca 11 IR triplet and a negligible impact on the atomic polarization 
of the upper levels of all the above-mentioned spectral lines. 

We have also shown that the collisional transfer rates have an unnoticeable impact on the 
emergent fractional linear polarization for most of the lines of multiplet 42 of Ti i, except per- 
haps for the /i4536 line of Ti i for which it is easier to detect a 25% reduction in the emergent 
polarization amplitude. This is important because such a line is insensitive to the presence of 
magnetic fields, which implies that we can use its scattering polarization pattern for investigating 
the reliability of 3D (magneto)-hydrodynamical models of the solar photosphere. 

CoUisional transfer rates might be truly important for other atomic or molecular systems 
with closer 7-levels (e.g., hyperfined structured multiplets and/or molecules). In fact, there are 
both theoretical and empirical indications that they could be significant for molecules like SiO 
(Derouich 2006) and MgH (Asensio Ramos & TrujiUo Bueno 2005; Trujillo Bueno et al. 2006). 
Therefore, future theoretical and experimental research in this direction will be worthwhile. 
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Appendix A: Multiplet 42 of neutral Titanium 

In what follows only the elastic and inelastic coUisional rates (in s"') are given as a function of 
neutral hydrogen number density uh (in cm"^) and kinetic temperature T (in Kelvins). However, 
it is straightforward to retrieve the values of the superelastic coUisional rates C^^\aiJi, a^Ju) by 
applying the relation: 

fin 23 1 + 1 Ea J ~ ^aiji (K) 

Q '{aiJuOuJu) = 7 exp — C) 'iauJu,aiJi), 

IJu + 1 KB-/ 

Eaj being the energy of the level {aJ) and kg is the Boltzmann constant. 

(K\ 

To simplify the notation we omit writting a in the coUisional rate symbol. Note that C) ( J, //) 
denotes the multipolar component of the inelastic coUisional rates that induce transitions between 
a lower level \JiMi > and the level \ JM > (both belonging to the same spectral term). Note that 
the largest population transfer rate is obtained for the lowest // value and that the semi-classical 
theory gives zero coUisional transfer rates for /S.J — J — 7/>2. With obvious notation, these are 
the colUsional rates for the Ti i levels of Fig. 1 : 

AA. Elastic coUisional rates: lower levels 



£)(2)(1) 


= 8.69 X 10 


"^^5000^ 


s~ 




= 8.35 X 10" 


"^^5000^ 


s~ 




= 8.19 X 10" 


-%h( ^ )0-4i 
5000 


s~ 


£)(2)(4) 


= 7.35 X 10 


"^^000^ 


s~' 




= 7.46 X 10" 


"^^5000^ 


s~ 



A.2. Population and alignment transfer rates: lower levels 



Cf\2, 1) 


= 11.04x10" 


"^^5000^ ' 


Cf\2, 1) 


= 3.22x10"' 


-^5o'oO>"-^^^- 


Cf\3, 1) 


= 5.31 X 10"'' 




Cf\3, 1) 


= 4.37 X 10"' 




cf \3,2) 


= 3.05x10"' 




cf \3,2) 


= 2.01 X 10"' 




Cf\4,2) 


= 2.91 X 10"' 


-^5o'oO^°"^- 
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Cf '(4,2) = 2.77 X 10"'^ tini 
(4,3) = 3.02 X 10"'' HHi 
(4,3) = 2.99 X 10"^ HHi 



5000 



Cf (5,3) = 2.95 X 10"^' nni 

(5,3) = 2.88 X 10"^ HHi 

Cf'(5,4) = 3.33x10-%h( 

Cf\5,4) = 3.12 X 10"'* HHi 



-9 



5000 



5000 



5000 



5000 



)0.42 ^- 



^0.43 ^- 



)0.42 



)0.44 ^- 



xO.42 



5000' 



5000 



A.3. Elastic collisional rates: upper levels 



D^^\l) = 1.21 X 10"'^ riHi 



5000' 



D<2)(2) = 1.00 X 10-^ "«^5^^"*' 

D<2)(3) = 0.81 X 10-9 "H^s^)"'^' 

£)(2)(4) = 0.93 X 10'' nni^f'' s 

D^'\5) = l.llxlO-%H(3^)°-^s- 



A.4. Population and alignment transfer rates: upper levels 



C 

d 

c\ 
c\ 

'-I 

c[ 
c\ 
C 
c 



''(2,1) 
'(2,1) 
''(3, 1) 
'(3,1) 
''(3,2) 
'(3,2) 
''(4,2) 
'(4,2) 
''(4,3) 
'(4,3) 
''(5,3) 
'(5,3) 
''(5,4) 



8.54 X 10" 
0.54 X 10" 
1.51 X 10" 
1.37X 10" 
1.90 X 10" 
0.65 X 10" 
1.16X 10" 
1.15x10" 
1.72 X 10" 
0.99 X 10" 
0.93 X 10" 
0.94 X 10" 
1.40 X 10" 
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Cf\5A) - 1.02xlO-%«(^rS-i 



Appendix B: Ionized calcium 

With obvious notation, these are the collisional rates for the Ca n levels of Fig. 5: 
6. 1. Elastic collisional rates 



D'-\Dy2) = 1.99 X 10-' "hC^)""" (B.l) 
D^^\D5/2) = 2.11 X 10 " nni-^f s'^ 



B.2. Population and alignment transfer rates 



cf\P3/2,Pm) = 5-70X io-%«(^r' 

cf\Ds,2,Dy2) = 2.21 X 10-' nH(^r' s-1 

Cf\D5/2,Dy2) = 1.06 X 10-' 
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